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1. Introduction 

Over the last year, the connection between anti-de Sitter (AdS) spaces and conformal field 
theories (CFTs) [1-3], has provided a method to study strongly coupled field theories. In 
gauge theories, a natural observable is the Wilson loop. In the gravity description, these 
are related to string worldsheets ending on the boundary of AdS. In this paper, we show 
how to calculate the operator product expansion (OPE) of the Wilson loop in which we can 
approximate a Wilson loop by local operators when it is small. Since we are dealing with 
a conformal field theory, the Wilson loop should be small compared to the distances which 
separate it from other Wilson loops or other operators. 

Our strategy for computing the OPE of a small Wilson loop in the 4D SYM theory is 
as follows. At large N and g Y M-^\ the bulk IIB string theory is in the small string tension 
and small curvature limit, so that classical string theory is a good approximation. In this 
context, the loops are represented in the bulk by classical (minimal area) worldsheets which 
end on the AdS boundary [4,5]. For example, two circular Wilson loops of radius a, which 
are separated by a distance L, correspond in the bulk to a worldsheet with these two loops 
as its boundary. When the ratio of the size of the loops to their separation is very small, 
a/L <^ 1, the worldsheet degenerates into two hemispheres connected by a very thin tube [6]. 
This degenerate worldsheet represents the exchange of light degrees of freedom in the bulk 
between two otherwise unaffected minimal surfaces, each with a circle as its boundary. It 
is then straightforward, albeit slightly tedious, to extract the OPE by properly identifying 
the light states being exchanged and their coupling to the strings. For Wilson surfaces, the 
strategy is the same, though one now considers membrane worldvolumes in M-Theory rather 
than strings. Here we find a new divergent term in the calculation which is proportional to 
the rigid string action. 

The paper is organized as follows. In section 2, we examine the OPE for a circular Wilson 
loop in some detail. We discuss what operators are allowed to appear in the OPE and list the 
operators of low conformal weight explicitly. We compare this with the leading terms in the 
perturbative expansion of the Wilson loop. We then discuss two approaches to the problem 
of calculating the coefficients of the OPE at large coupling. The first method consists in 
calculating the correlation function between a Wilson loop and the various operators of the 
theory, (W(C)Oi). The second method is to compute the correlator between two Wilson loops 
and to identify the contributions at each order in the size-to-separation, a/L, expansion. 

In section 3, we find the minimal area string worldsheet that describes a circular Wilson 
loop in the fundamental representation. In section 4, we identify the scalar modes that 
contribute to the supergravity interaction between two Wilson loops and consider their cou- 
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pling to the string worldsheets. We then compute the Wilson loop correlators and extract 
the OPE coefficients. In addition we consider the potential between two rectangular Wilson 
loops, which is a straightforward application of the same techniques used to compute cor- 
relators of loops. In section 4.4, we outline some qualitative features of the computation of 
exchange of tensor modes between the worldsheets. 

In section 5, we consider Wilson surfaces in the (0,2), d = 6 superconformal theory. We 
find the minimal area membrane worldsheet solution describing spherical Wilson surfaces 
and the supergravity modes of AdS 7 x S* 4 which contribute to their correlation. We find 
that there is a UV logarithmically divergent contribution to the area of the surface that is 
proportional to the action of a rigid string embedded in the 5-brane worldvolume. 

Several details about the bulk-to-bulk and bulk-to-boundary Green's functions that we 
will need are presented in an appendix. 

2. The Operator Product Expansion of the Wilson Loop 

In [4,5], a prescription was given to compute the effective quark-antiquark potential in the 
large N strong coupling limit of maximally supersymmetric Yang-Mills (conformal) field 
theory in four dimensions. This potential could be obtained by computing the expectation 
value of the Wilson loop, 



where I is the total length of the loop and $ is a UV regulator [7]. In this paper, we consider 
the modified Wilson loop operator given by 



where A^a) are the four- dimensional gauge fields, the X 7 (<t), I — 1, 6 are the six scalar 
fields of Af = 4 SYM, and # 7 (s) is a point on the five-sphere (so 9 2 = 1). As argued in [4,5] 
this is the operator which leads to a simple calculation in supergravity. The reason is that, 
if the gauge symmetry is broken to U(N — 1) x U(l) through a Higgs VEV, the massive 
W-bosons can be interpreted as strings stretched between the horizon and a D3-brane in 
AdS. These W-bosons carry, in addition to the charge under the gauge fields, a "scalar" 
charge under the scalar X I Q 1 ', where 1 is the 5*0(6) orientation of the Higgs VEV, so that 
we get the second term in (2.2). 

We expect that there exists an operator product expansion for the Wilson loop when it 



(W(C)) ~ lim e~^- m \ 



(2.1) 




(2.2) 
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is probed from distances large compared to its characteristic size a, 



W(C) = (W{C)) 



E(n) A (n K n (n) 
c\ 'a ' O ' 



(2.3) 



where the are a set of operators with conformal weights A\ n \ In this notation, we let 
Of^ denote the i th primary field, while the for n > are its conformal descendants. 
For the circular Wilson loop solution, the expectation value of all operators other than the 
identity vanishes, so that the coefficient of the identity is the expectation value of the loop. 

The problem is to explicitly calculate the coefficients cf 1 ^ that appear in the Wilson loop 
OPE. In the field theory, this can be done perturbatively at weak coupling, but, barring 
the existence of nonrenormalization theorems, the result cannot, in general, be reliably 
extrapolated to strong coupling. In fact we will see that the coefficients are different. 

Local gauge invariant operators are given by traces (over the gauge group) of polynomials 
of the scalars X 1 , the fermions ip a , the Yang-Mills field strength F^ u , and their covariant 
derivatives. Since the operators appearing in the Wilson loop must have the same symmetry 
properties as the Wilson loop itself, the operators should be bosonic and gauge invariant. 
We consider the case where # 7 (s) = 1 is a constant. This breaks the SO (6) R-symmetry of 
the superconformal field theory to 5*0(5), so we expect that the Of 1 ^ are SO (5) invariant. 
Therefore we will obtain operators which are in irreducible representations of SO (6) that 
contain singlets under the SO (5) maximal subgroup. 

We perform the analysis at each conformal dimension. 

• A = 0. The only operator of dimension zero is proportional to the identity. 

• A = 1. The only elementary fields of dimension one are the scalars X 1 . A gauge 
invariant operator would be tr(X 7 ), but taking SU(N) as the gauge group [3], this 
trace vanishes. 

• A = 2. At dimension two, the only trace over the gauge group that is non-trivial is that 
of the scalar bilinear tr(X 7 X J ), which splits in two irreducible representations of 5*0(6). 
One is the singlet tr(X / ) 2 which is expected to have a large anomalous dimension [3], 
so it must be dropped from the expansion at low orders in the supergravity regime. 
The other is the symmetric traceless tensor, Cf J tr(X I X J ), which is the 20 of 50(6). 
Its conformal dimension is independent of the coupling [3]. Since the 20 — > 1 © 5 © 14 
under 50(5), we find that the operator in the 50(5) singlet, i.e., the projection of 
Cf J tr(X I X J ) will appear in the OPE. 
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• A = 3. At dimension three, we must consider the scalar operator tr(X I X J X K ). 
Once again, only the symmetric traceless part, Cf JK ti(X I X J X K ) (in the 50), has a 
protected conformal dimension. All other components ought to have large anomalous 
scaling dimensions. 

We must also consider the scalar bilinear in the fermions, tr (ip^ipBa), where A and B 
are spinor indices for the R-symmetry group. This operator is in the 10 — > 10 and 
does not contain any SO (5) singlets, so it will not contribute. 

We can also have tr(X / F M ^), which transforms as a two form under Lorentz transfor- 
mations and decomposes as the 6 — > 1 © 5 under 5*0(5) and should contribute to the 
OPE. For a circular Wilson loop, the allowed components depend on the orientation 
of the loop. 

The final primaries at this order are the R-symmetry currents, jj! J \ which are in the 
adjoint of 50(6), which is the antisymmetric tensor, 15. Under 50(5), 15 — > 5 © 10, 
so there is no 50 (5)-invariant component. Therefore this operator does not contribute 
to the OPE. 

Finally, we can also have Cf J d i tr(X I X J ), which is a superconformal descendant of 
Cf J tr(X I X J ). In the particular case of a circular Wilson loop OPE, it is forbidden by 
rotational invariance. 

• A = 4. At dimension four, there are various chiral primaries. The operators which 
contain 50(5) singlets and should appear in the OPE are the symmetric traceless 
rank 4 tensor Cfj KL tr(X 2 ' X J ' X K X L ), the field strength operator Cf J ti(X I X J F IJtv ) , the 
energy-momentum tensor, and the Lagrangian. Additionally, one can have descendents 
of the A = 2, 3 operators which already appeared above, as well as two trace operators 
like Cf J tr(X I X J )C^ L tr(X K X L ). In the 't Hooft limit, we expect to find only single 
trace operators. 



Summing up our results, we arrive at the following expression for the circular Wilson 
loop OPE 



where a^ u is a unit two-form which denotes the orientation of the Wilson loop, Af k is a 
constant that ensures that operators are "unit" normalized, in a sense indicated below and 
Cj^.j is a basis of symmetric traceless tensors such that the spherical harmonics are Y A (9) = 



W(C) 

WW) 



1 + cf a 2 Y%\6)N 2 CtM {X f X J ) + cf a 3 Yf (6)M 3 Cf JK tr {X 1 X J X K ) 
+ cf ) a 3 a^tr(e I X I F^) + -- - , 



(2.4) 
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Of r ..j 9 11 ■ ■ -9 Ik , with the index A running over all the spherical harmonics of given SO (6) 
Casimir, see [8] for conventions. 

2.1. Perturbative Calculation of the OPE Coefficients 

For small A = gy M N , we can perturbatively expand the Wilson loop (2.2) to find an expres- 
sion for the OPE (2.4). We find the first few terms to be 

- ^^-tr(^) - ^(-^) 2 + " " " • 
Some of these operators will get high conformal dimensions in the strong coupling limit and 
so we know that they will not appear as leading terms in the expansion. There are, however, 
operators whose dimensions are protected, such as the symmetric traceless combinations 
Cf x T tr(X Il ...X Ik ). The operator product coefficients will depend on the normalization of 
the operators. We will choose them to be "unit" normalized, in the sense that (0(x)0(y)) = 
\x — y\~ 2A . The operators in (2.5) are not normalized. In order to normalize them, we have 
to compute their two-point functions. These were calculated in [8] and using those results, 
we find that the operator product coefficients for the highest weight chiral primaries have 
the behavior 

A A/2 

ca ~ — (2.6) 

where A = gy M N . For the other protected operators, we find a similar \ P /N dependence, 
where the exponent p is related to the conformal weight of the operator. We will see that 
for large A, the dependence on N will be the same, but that the A dependence is different. 
Of course the dependence on N can be understood in a simple fashion by using large N 
counting arguments. 

2.2. Supergravity Calculation of the OPE Coefficients 

We now describe how to calculate the coefficients in the supergravity description. There are 
two ways to determine the coefficients of the operator product expansion (2.4). The first 
and most straightforward method is to compute the correlator of the Wilson loop with each 
operator that is expected to appear in the loop. This correlator gets contributions only from 
the given conformal primary and its descendents, 

(W(C)Of ] ) _ (o) a A<i ' r M„A^/ r9 (mWon (2 7) 
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Here we have isolated the contribution from the descendents and their mixings with the 
primaries in the second term. 

In the supergravity description, the Wilson loop will be related to a string worldsheet 
ending on the boundary of AdS. The correlation function (2.7) can be calculated by treating 
the string as an external source for the fields propagating in anti-de Sitter spacetime and 
then computing the string effective action for the emission of supergravity states onto the 
point on the boundary where the operator is inserted. See Figure 1. 



xi 



x 2 




z=0 



Figure 1: The emission of a particle from a circular Wilson loop of size a onto the AdS 
boundary at a distance r from the loop. 

Another approach to the problem of calculating the OPE coefficients is to compute the 
correlator of a pair of Wilson loops that are separated by a distance which is large compared 
to their size. In the conformal field theory, the correlator can be calculated from the operator 
product expansion for the two Wilson loops 

(W(C, L)W(C, 0)) m a(-).la(«) 



(W(C,L))(W(C,0)) 



£ c^cf a^'^r (0t\L)0f } (0)> 



;,j;m,n 



E(4 0) ) 



(o)\2 a 



2A 



(0) 



L 2A 



(o) 



(2- 



(m) (n) A (m) +A (n) 



{0 <r\L)o<r\o)). 



In the last line, the first term is due solely to the primary fields, while the second contains 
the contributions from descendents. 

In the supergravity approximation, these Wilson loop correlators can be calculated by 
computing the amplitude for the exchange of light states between the two string worldsheets 
which have the Wilson loops as their boundaries [6], as represented in Figure 2. We will 
actually calculate the OPE coefficients in this fashion, since it is slightly simpler. 
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z=0 



Figure 2: The correlation between two circular Wilson loops. 

In the next section, we will examine the details of the circular Wilson loop solution. With 
the necessary information in hand, we will then return to the computation of correlation 
functions of Wilson loops. 

3. Circular Wilson Loops and AdS Supergravity 

According to [4], in order to compute the expectation value of a circular Wilson loop in the 
large g s N limit, we should find the minimal area string worldsheet ending on a circle at 
the boundary of anti-de Sitter space. We choose the scalar charge of the Wilson loop to be 
constant, so that the string worldsheet lives at a single point on the 5-sphere. This implies 
that # 7 (s) = 6 1 in (2.2) is a constant. We could find the classical worldsheet by solving the 
Euler-Lagrange equations coming from the Nambu-Goto action in this background, however, 
in this case there is a simpler way to find the worldsheet. 

We note that the Euclidean conformal group in 4 dimensions, 5*0(1,5), has elements 
which map straight lines into circles, namely the special conformal transformations, 

a _ x l + dx 2 

l + 2c-x + dx^ 1 } 

where d is a vector in IR 4 . We take the AdS metric 

ds 2 = -4 (dz 2 + dx l dx l ) , (3.2) 
with the boundary at z — 0. The special conformal transformation (3.1) corresponds to the 
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AdS reparameterization 



H x % + c l (x 2 + z 2 ) 

X — 



1 + 2c-x + c 2 (x 2 + Z 2 ) ,ooA 
z [6 - 6) 



z = 



1 + 2c • x + c 2 (x 2 + z 2 ) ' 



Starting with a line on the boundary and a minimal area surface in AdS which ends on 
that line (which just extends along the line and the z direction), we can apply the above 
conformal transformation and map it to a circle in such a way that on the boundary we have 
a circle of radius a. Then the surface in AdS is given by 1 . 

X H = Va 2 - z 2 (ej cosip + e* 2 sinip) , < z < a, < ip < 2vr. (3.4) 

where e±, €2 are two orthonormal vectors on the boundary. We see that this surface ends on 
the boundary on a circle of radius a and it closes off at z = a. It is useful to compute the 
area element, 

dA = dzdip^Jdet (g^(x)d a x' 1 dpx u ) = a< ^ z ^ _ 



As in [4], we find that there is a divergence in the action, 



S e = 

2na 



1 r , f 2n , , az 1 (a \ 



In terms of the theory on the boundary, this divergence corresponds to the UV divergence 
in the Coulombic self-energy of a point charge. In the bulk theory, this divergence is due to 
the contribution of an infinitely long straight string ending on the circle. After subtracting 
the divergent term, we find that S = —1/a', which is independent of the radius a of the 
loop, as required by conformal invariance. We are choosing units in which the radius of AdS 
is equal to one, so that a' — 1/ \/Aixg s N . 



4. Contributions from the Lightest Scalars 

We will be primarily interested in the contributions from the lightest scalars, whose exchange 

will dominate the long distance interactions. These light states correspond to the operators 

of lowest dimension in the OPE for the Wilson loop. The relevant modes may be determined 

from the KK mass spectrum listed in Figure 2 or Table III of [10]. 

1 While this paper was being written, we learned that D. Gross has independently obtained this circular 
Wilson loop solution [9]. 



8 



4.1. The Dilaton 



We will start calculating the coefficient of the OPE in front of the operator associated to the 
dilaton. We start with this field because the calculation of the precise coefficient is simpler. 
We will consider other cases later. The dilaton can be expanded in Kaluza-Klein harmonics 

as 



= £^y<*>(0). 



(4.1) 



The action for the dilaton is 



/^£i^[(V0*) 2 + m^ 2 



(4.2) 



where, in units where i?Ads — 1) 



k 2 fc - 1 7r 2 (A; + l)(A; + 2)' 



(4.3) 



and we are normalizing the spherical harmonics as in [8]. 

The coupling of the dilaton to the string worldsheet can be found remembering that the 
supergravity calculations were done using the Einstein metric, while the string worldsheet 
couples to the string metric, G str ing = e^ 2 g Einstein- So the coupling of the string worldsheet 
to the dilaton is given by 



'dii 



(4.4) 



where 7^,3 is the induced metric on the worldsheet and is the corresponding worldsheet 
curvature. In the case that the radius of AdS is large, the term involving the worldsheet 
curvature will be subleading compared to the first term. Therefore we should neglect the 
curvature contribution in the amplitudes we consider. 

So now we want to calculate the contribution of the dilaton to the expectation value of 
the product of Wilson loops. This contribution will be given by 



{WW) 



dil 



(wy 



Exp 



k,A 



(4.5) 



where the sum over A indicates the sum over all spherical harmonics of total angular mo- 
mentum J 2 = k(k + 4) 2 . This formula summarizes the effects of dilaton exchanges between 

2 This sum gives a result that is independent of 8, but we leave the result in this form to read off the 
contribution to the OPE of each Kaluza-Klein mode. 
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the two worldsheets. Of course, the dominant term for large N is the one dilaton exchange 
arising from the first-order expansion of the exponential. Notice that a, a' indicate points 
on the two separated worldsheets. If the separation between the worldsheets is very large, 
then W is very small and we can approximate the Green's function (A. 9) by 

G t (WU ~ ~ (4.6) 

where L is the (large) distance between the two loops, A = 4 + k and a = (A — l)/2ir 2 . Us- 
ing (3.5) and (4.6) we find that the integrals over the worldsheet reduce to simple expressions 
involving J dz z A ~ 2 which are always convergent at z — 0. The final result is that 

(WW) dil ^ 2 k ~ 2 7r(k +l)(k + 2) g s ,. Q . 2 / a\ 2A 



E V + 3 J l^' 2 (l) < 4 ' 7 > 



so that 



Cd«,A = 2 y jj- (4.8) 

Notice that we are unit-normalizing the dilaton operator. 

Of course we could have calculated this last result directly by computing the one point 
function of the operator associated to the dilaton in the presence of a Wilson loop. In fact, 
as shown in [11], the one point function is proportional to the value of the dilaton near the 
boundary. This value is non-zero because the string worldsheet acts a source for the dilaton. 
We can actually calculate the correlation function of the operator with the Wilson loop for 
any position, not just for large distances. We get 

{WO\x)) _ n v _ _ fc , . _ , ^ . 1 r 



c dUA Y {k) (9) 



a 



2A 



(4.9) 



[(y 2 + r 2 - a 2 ) 2 + 4aV] A/2 ' 



where r is the polar coordinate on the plane defined by the loop and y is the distance on 
the plane orthogonal to the loop. We see that when the operator approaches the loop we 
have a singularity of the form l/d A where d is the distance from the loop. In fact, the 
result (4.9) could be derived in a simple way by first computing (WO) for a straight line and 
then applying the conformal transformation mapping the line to a circle. This result (4.9) 
includes all the information about the conformal descendents of the operator O appearing 
in the OPE. 
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4.2. "Tachyonic" Scalars 

The leading term in the OPE expansion of the Wilson loop comes from an operator of 
dimension two. This operator comes from a field with m 2 < in the supergravity theory. 
Scalars arise from several supergravity fields. There are scalar KK modes of the metric over 
the 5- sphere, (in the notation of [10]) 

K(x,6)=J2* k (x)Y ik) (6), (4.10) 

k 

as well as the scalar KK modes of the antisymmetric 4-form, 

Y,b k e a ^D f Y ik \B). (4.11) 

k 

One linear combination of these has tachyons in its spectrum. Another scalar comes from 
the trace of the metric on the AdS component, M. We can algebraically express these in 
terms of the ir k , as 

16 ( 412 ) 

k 

so they contribute to the states of this family. 

From the field equations, given as equation (2.33) of [10], one sees that the modes ir k and 
b k mix. The mixing angles, as well as the normalized action for the mass eigenstates has 
been conveniently presented by Lee et. al. [8]. They found that the mass eigenstates were 

s k = 1 [n k - W(k + A)b k ] , = k(k - 4), k > 2, 

2U(/£ + 2j (4.13) 
t k = 2Q{k + 2) + Wkb k ] , Ml = (k + 4)(fc + 8), k > 0. 

The lightest states correspond to the lowest modes of s k , so we will focus on these modes. 
The action for the s k was found to be [8] 

S= I <£xY,B k l -[{V,s k Y + k(k - 4)(s k ) 2 } , (4.14) 



where 



2*- k N*k(k - 1) 

Bk= **(k+l)* ' (415) 
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and our normalizations are as in [8]. 

In order to compute the coupling of s to the string worldsheet, we need to find all the 
supergravity fields that are excited when s is nonzero and all the other "diagonal" modes are 
set to zero. From the equations that we saw above we find a contribution to and there is 
also a contribution to 

H w = kTi D <* D ") 8 > (416) 

where the parenthesis indicate the symmetric traceless combination. 

Now we should find how s couples to a string worldsheet. These couplings will involve 
terms with derivatives. In the calculation we are interested in these derivatives will act 
on the Green's function of the field s. Moreover, we will be interested in extracting the 
leading piece in L, where L is the separation between two Wilson loops or a Wilson loop 
and an operator. In that regime we will be able to approximate the Green's function by an 
expression like z A /[(x — x') 2 + z 2 ] A , so that the only derivatives which will not produce a 
subleading term in 1/L will be those acting on the numerator of this expression. This implies 
that in calculating the coupling to s from (4.16) we will be able to replace derivatives by 
factors of A, etc. The string worldsheet will couple to the various components of the metric 
on AdS$. Couplings through a a ^s will be small if R is large, since they will involve couplings 
to the worldsheet fermions. Finally we get the coupling to the worldsheet as 



2na 

By using the same method as we used for the dilaton we obtain 
(WW), « ( 2A; 



^)^W^£^i>W^. (4.18) 



{W) 2 , ( „ x , 

Similarly, we can calculate 



From these expressions we determine the OPE coefficients 

c. >A = 2*/ 2 -VA^. (4.20) 

This equation should be compared to the weak coupling result (2.6). As expected the iV 
dependence is the same but the powers of Qym^ are different. This is no contradiction, since 
the two calculation have different regimes of validity. 
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4.3. The Potential Between Two Rectangular Wilson Loops 

The tools that we have collected to compute correlation functions due to exchange of scalar 
supergravity modes between string worldsheets are also applicable to the study of the po- 
tential between rectangular Wilson loops. For a pair of rectangular Wilson loops, each of 
size a, which are separated by a distance L ^> a, as depicted in Figure 3, the potential takes 
the form [12] 




a 



Figure 3: The interaction between two rectangular Wilson loops of size a, separated by a 
distance L ^> a. 



v = _ hm Mg. L)W{C, °)) - Mg, L))(W(C, 0)) 

,2n-2 (4.21) 



W- 

Vn T,2n-1 ■ 



L 2r 

n 

Perturbatively, we will find that the V n ~ X p /N 2 . At large A, the iV dependence is the same, 
but the A dependence will be different. 

We will use the worldsheet solutions of [4] to compute the effective potential. We find 
that, for large separations, the asymptotic behavior of the scalar Green's function (A. 9) is 

WJ-^^y.^V--- (4-22) 



The surface was given by 



•'•(■-)= '-wo* I dy y 4 , (4.23) 

z/Zmax V 1 — U 
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where z max is determined by the condition that \x(0) — x(z n 
ar(l/4) 2 /(27r) 3/2 . The area element is given by 



dA = dtdz 



The coupling to the field s is given by 



1 



2na' 



dA- 



'-2ks)z 2 



= a/2, so that z max = 
(4.24) 

(4.25) 



where we used the same method as above to calculate the coupling involving terms with 
derivatives of s. So the final expression for the potential between two Wilson loops is 



= - 



T 

y- a T(k- 1/2)0F 



-I 2 



2na' 



(4.26) 



l\4fc-4 



322fe+5 7r 3fe-7/2 



fc(fc + 1) 2 (A; + 2)(fc + 3)r(fc - |)r(^) 4 ^ s iV q 2fc ~ 2 

iv 2 L 2k 1 ' 



; - . r(A;)r(^i) 2 

We have done the sum J2a y a( ) 2 = (k + 3)(k + 2)/(3-2 k+1 ). We have included contributions 
to the potential coming from all of the Kaluza-Klein modes of the field s. The ones with 
k = 2 will give the leading contribution. Of course other fields will also make contributions 
to the potential which are comparable to the terms above with k > 3. For example, we can 
calculate the contributions to the potential from the dilaton 

2 



k,m k 



a r(A-l/2)0F 



^ t(A)l 2 a-i 



A-l- 



2ixa' 



4r 



A+l* 



r 



l\4A-4 



(A - 1) 2 (A - 2) 2 (A - 3)r(A - \)T{^Y 



322A+97J-3A-7/2 



r(A)r(^i) 2 



(4.27) 



g s N a 



2A-2 



N 2 L 2A_1 ' 

where now A = 4 + k. As g s ~ A/iV, these results agree with the behavior illustrated 
in (4.21). 



4.4. Contributions from Vectors and Tensors 

Let us finally note that the operator product coefficients involving other operators can be 
computed in a similar way. They will involve the contribution of various supergravity fields 
to the correlator between two different Wilson loops. In particular, the term tr(X I F IM/ ) term 
in the OPE (2.4) corresponds to the lowest mode of an antisymmetric tensor, B^, on AdS 5 . 
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5. The Spherical Wilson Surface and AdSy 



In [4], it was shown that one could use the AdS description of the large N limit of the (0, 2) 
superconformal field theory in six dimensions to compute Wilson surface observables [13], 
even though an explicit formulation of the field theory does not exist, so that there is no 
formula analogous to (2.2). Let us consider a spherical Wilson surface. We take the scalar 
charge of the surface to be constant (a point on S 4 ). In the gravity picture, the Wilson 
surface should be the boundary of a minimal area membrane worldvolume in AdSy x S 4 . 
One can either solve the equations of motion directly to obtain the minimal worldvolume, or, 
by analogy with the discussion of section 3, one can note that a flat plane in the boundary of 
AdSy can be conformally mapped to a sphere. This flat plane is the boundary of an infinite 
membrane that is stretched between the AdS boundary and z = oo. The conformal mapping 
maps the worldvolume into a 3-hemisphere whose boundary is a 2-sphere that corresponds 
to the CFT Wilson surface. A convenient parameterization of the solution is given in terms 
of the Poincare coordinates as 



'a 2 


-z 2 


'a? 


-z 2 



x 2 = V a 2 — z 2 sin 9 cos ip (5.1) 



x 3 = V 'a 2 — z 2 sin# sin^, 

where < z < a, < 9 < ir, and < ip < 2n. Now we take the radius of AdS to be equal 
to one, then the radius of Rs 4 = 1/2, l v = 1/(8tcN) 1 ^ 3 , and the tension of the two-brane is 
= 1/(2tt) 2 /3 = 2N/n. 
We then find that the volume of the membrane is divergent 



a dzay/a 2 — z 2 



S = T® J dV = T^47fJ^ ^ 



7° 

( 5 - 2 ) 



2a 2 , 2a 
V— -2 In— -l + 0(e 



e 



We can make several observations regarding this expression for the action. First, we see that 
the action scales as N, in agreement with the scaling found for the "rectangular" solution 
of [4] . As indicated in [7] , e should be thought of as a UV cutoff. So we see that we have two 
divergent terms. The quadratic divergence is proportional to the area of the surface. This 
term was also present in the case of a rectangular Wilson surface [4] . In this case, we see that 
there is also a logarithmic divergence. The first question is what would this divergence be 
in a more general case? It can be seen, by analyzing the equations of motion of the theory, 
that for a generic two-dimensional surface S 2 this logarithmic divergence is proportional to 
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the "rigid string" action [14] 

Srigid = [ dVV7(V 2 ^) 2 , (5.3) 

where 7 is the induced metric on the Wilson surface and X 1 are the coordinates on R 6 
describing the surface. Notice that S 2 is a surface in the boundary six-dimensional field 
theory. As emphasized in [14], this action is invariant under scale transformations in the 
target space, X 1 — > XX\ which is consistent with what we expect in a conformal theory. 
Actually, it is possible to prove that the action is also invariant under special conformal 
transformations, so that it is invariant under the full conformal group 3 . One implication of 
this logarithmic term is that the expectation value of the Wilson surface is not well defined, 
since we can add any constant to the logarithmic subtraction. Furthermore, it seems to 
indicate that the expectation value of a Wilson surface is scale dependent. 

It seems natural to speculate that tensionless strings in this six-dimensional field theory 
are governed by some supersymmetric form of the the action (5.3). 

Despite the fact that the expectation value of the spherical Wilson surface is not well 
defined, the connected correlation functions of Wilson surfaces do not receive extra divergent 
contributions. These correlators can be calculated in a completely analogous fashion to the 
Wilson loops in section 4. One considers a Wilson surface whose characteristic size is much 
smaller than its distance from any probe in the theory. Then one identifies the operators 
that are allowed to appear in the OPE and computes the necessary correlation functions to 
extract the OPE coefficients. The details will appear in [15]. 

6. Conclusions 

In this paper, we have made use of the the connection between Anti-de Sitter spaces and 

conformal field theories in the large N and gy M iV limit to compute the operator product 

expansion of "small" Wilson loops in H — 4, d — 4 super Yang-Mills theory and Wilson 

surfaces in the (0, 2) superconformal theory in six dimensions. 

By determining what supergravity states couple to the worldsheet describing the Wilson 

loop, we found the set of operators that are allowed to appear in the OPE. By computing the 

amplitudes for exchange of supergravity modes between a Wilson loop and the boundary and 

between two Wilson loops, we were able to compute the correlation functions of a Wilson 

3 It is enough to prove that the action is invariant under inversions X' 1 — ► X' 1 / X 2 . This can be shown 
using identities like V pX l \7 2 X l = 0, which use the fact that -f a p — d a X l dpX l is the induced metric. 
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loop with a CFT operator and with another loop. From these expressions, we were able to 
deduce the coefficients that appear in the OPE for the states that we considered. 

We also investigated Wilson surfaces in the (0,2) six-dimensional field theory and we 
found that there is a UV logarithmically divergent contribution to its expectation value 
proportional to the action of a rigid string embedded in the 5-brane worldvolume. 
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A. Green's Functions 

Scalar Green's functions on anti-de Sitter space have been discussed in a large number of 
papers [3, 2, 16-19] in connection with the types of correlation functions we are looking to 
compute. We will repeat some of these calculations in this appendix. 

Consider the action of a real scalar field on Euclidean anti- de Sitter spacetime, of 
radius one, with source J, 



d d+1 Xy/gB 

AdS 



\m 2 + ^ 2 -<i>j 



(A.l) 



where B is some constant. The equation of motion is 

B{-Vl + m 2 )<P = J, (A.2) 

where V 2 , = -^d^y/gg^d^) is the Laplacian. Imposing the boundary condition (f)\dM = 
yields a unique solution for 0, as long as the operator — V 2 + m 2 , is positive definite. This 
is the case for all m 2 > -d 2 /A [20]. 

Solutions for which minimize the action (A.l) are given by the integral equation 

<f>(x)= [ d^x'^gJ^Gix^x^Jix'), (A.3) 

J M 
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where the kernel G(x,x') is the covariant Green's function for the equation of motion (A. 2), 
satisfying 

B(-V 2 X + m 2 )G(x, x') = * 5 {d+1 \x - x'). (A.4) 

By evaluating the action at a solution of (A. 3), one finds that 

S = -\f d d+l x [ d d+1 x' yfij{x)J(x)G(x,x!) y/ffijjtf). (A.5) 

2 J M JM' 

If we consider the specific case of two separated sources, J±(x) and ^(x'), the amplitude for 
their interaction is then given by 

S=- I d d+1 x [ d d+1 x' ^{x)J 1 (x)G(x,x') y/gWMxf). (A.6) 

JM JM' 

A.l. Bulk-to— Bulk Scalar Green's Functions 

We will find it most convenient to compute the Green's function in the upper-half space 
representation of anti-de Sitter spacetime, with metric 



ds 2 = 

z 



^(dz^Y.dx^j. (A.7) 



As the scalar Green's function can only depend on the distance d(x, x') between the sources, 
in this metric it will be a function of 

W = - = - - (A8) 

{z-z>Y + Yi=i\ x i-<\ 2 2coshd(x,x')-l' 

Note that W is singular precisely at (z,x) = (z',x'), which is the location of the singularity 
in the Green's function (A.4). It is easy to show that the solution for (A.4) which goes to 
zero at the boundary is given in terms of the hypergeometric function 



G(W) = ^W A 2 Fi(A, A + 2A - d + 1; -AW). (A.9) 



where A = d/2 + ^m? + d 2 /4 will be the conformal weight of the associated operator and 
«o is 

r(A) 

i + (A. 10) 



Otn= — „ 

2vr fr(A-f + i) 



a = - 9 for d = A 
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A. 2. Computation of correlation functions 

Let us consider the field in the bulk (A. 3) produced by a source on the boundary. To be 
precise, we take J{x') to have support very close to the boundary, 

supp(J«)) = {z'\0 < z' < e}, (A.ll) 
for infinitesimal e. Then, using the bulk-to-bulk Green's function (A. 9), we can write 

<f>(x)= [ d d x' I dz' y/^)G(x,x')J(x'). (A. 12) 

JdM JO 

In the region over which we integrate z', W ~ 0, so that we can approximate 

G(W) ~ ^-W A , (A.13) 

B 



so that 

■ / d' 1 .,-' ( V" 



~ ^ J bm M ( , 2 + | g_ f? ) [dx' z^J(x'). (A.14) 



Now, to make contact with Witten's analysis in [3], we want to define a dimension d — A 
source 4>o(x'). 

o (f ) = c f dz' z'^^Jix'), (A.15) 
Jo 

where c is a numerical factor which we will determine from the 2-point function. Now, 
evaluating the action obtained from (A. 6), we find that 

S=^J d d xd d x' |f _ 1 f/|2A j dzz^- d - l J(x) j dz' z' A - d -\j(x) 

a ° f d d xd d x' ^°( f )^°( f/ ) ^ A ' 16 ^ 



\x — X' 



Be 2 

The two-point function is then 

5 2 S a 



5>o(x) £00(^0 Be 2 \x — x'\ 2A 



(A.17) 



Choosing the convention that the operator corresponding to this scalar is unit-normalized, 
(OO) = l/x 2A , we determine that c = ^a /B. 
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